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We discuss a kinetically constrained model in which real- valued local densities fluctuate in time, 
as introduced recently by Bertin, Bouchaud and Lequeux. We show how the phenomenology of this 
model can be reproduced by an effective theory of mobility excitations propagating in a disordered 
environment. Both excitations and probe particles have subdiffusive motion, characterised by dif- 
ferent exponents and operating on different time scales. We derive these exponents, showing that 
they depend continuously on one of the parameters of the model. 



I. INTRODUCTION 

There has been considerable recent interest in the hy- 
pothesis that glassy materials can be described by coarse- 
grained models with simple thermodynamic properties 
and non-trivial kinetic constraints P, H, H, H, H, S Hi 
These models capture the dynamically heterogeneous na- 
ture of glass-formers ^ : the implicit assumption is that 
microscopic details of the glass- former are important only 
insofar as they set the parameters of the coarse-grained 
dynamical theory. Some kinetically constrained models 
describe particles hopping on a lattic e [Tol [ill : in other 
cases, binary (Ising) spins are used [1, Il2l.ll3||. where the 
two states of the spin represent 'mobile' and 'immobile' 
regions of the liquid. 

In a recent paper, Bertin, Bouchaud and Lequeux 
(BBL) (5(1 discussed a kinetically constrained model in 
which molecular degrees of freedom are modeled by a 
real-valued local density, defined on a lattice. Loosely 
speaking, regions of high density correspond to immobile 
sites in the spin description of and regions of low 
density correspond to mobile sites. However, the contin- 
uous range of densities in the BBL model captures the 
fact that a glass-forming system has a variety of local 
packings, which may not permit a simple decomposition 
into mobile and immobile. The continuous range of den- 
sities leads to a continuous range of mobilities, resulting 
in a very broad distribution of relaxation times, charac- 
teristic of glassy behaviour. 

As discussed in Ref. [F], relaxation in the BBL model 
occurs by means of 'mobility excitations' that propagate 
subdiffusively across the system. Links between broadly 
distributed relaxation times and subdiffusive motion of 
particles are quite familiar in theories of glass-forming 
liquids [ij]. Here, we focus initially on the motion of 
mobility excitations, by coarse-graining the BBL model 
onto an effective theory for these excitations. This pro- 
cedure represents a very simple example of the coarse- 
graining of glassy materials that was proposed by Gar- 



rahan and Chandler 0. The resulting effective the- 
ory is a disordered generalisation of the one-spin facil- 
itated Fredrickson- Andersen (FA) model [l^- The pos- 
sibility of subdiffusive motion in this model is impor- 
tant when comparing model results with experiments on 
supercooled liquids: on approaching the experimental 
glass transition time scales increase dramatically, while 
the length scales associated with dynamical heterogeneity 
grow more slowly d, [l^ . In the dynamical facilita- 
tion picture 0, the motion of mobility excitations leads 
to a relation of the form t C^/f^ where z is a dynam- 
ical exponent, t the relaxation time scale, ^ the length 
scale associated with dynamical heterogeneity, and 11 a 
(possibly temperature dependent) constant [l^- Subd- 
iffusion of these excitations corresponds to an exponent 
z > 2, consistent with a time scale that increases much 
more quickly with the corresponding length scale than 
for ordinary diffusion. 

In this paper we focus throughout on the one- 
dimensional case, where subdiffusion effects are most pro- 
nounced [3] • Analysis of the disordered FA model leads us 
to two main results. Firstly, we are able to explain the 
scaling exponents observed in In particular, while 
the disorder in both the BBL and disordered FA mod- 
els is fluctuating, we explain why excitations propagate 
with the scaling laws expected for a particle moving in 
a quenched random environment. Secondly, we consider 
the motion of probe particles in the BBL and disordered 
FA models. These particles propagate subdiffusively, but 
with scaling laws that are different from those of the mo- 
bility excitations. 

The form of the paper is as follows. In section [Hi 
we define the BBL model and the disordered FA model. 
In section lllli we use four-point correlation functions 
[1, [H, [H, li^l to investigate the subdiffusive propaga- 
tion of mobility excitations, and we discuss the associated 
scaling exponents. In section HVl we consider the motion 
of probe particles in the BBL model, and show that this 
behaviour can also be reproduced in the disordered FA 
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model. 



II. MODELS 

A. BBL model 

The (one-dimensional) BBL model is defined Q for a 
chain of continuous densities {pi}, constrained to < 
Pi < 2. Dynamical moves involve rearrangement of the 
density between adjacent pairs of sites: 

p^, p^+i ^ p[, p[+-^, rate /jPp(p',,p^_^i|pj (1) 

where fi = 0(2 — — pi+i) is the facilitation function 
for bond i, between sites i and i + 1. Here, G(x) is 
the Heaviside step function, so a density rearragement 
between two sites can occur only if the total density on 
those sites is less than two. The distribution of densities 
after the rearrangement is 

Pp(p„ = A{p,p,+^Y-^5{R -p,- (2) 

where the delta function enforces volume conservation. 
The parameter p > Q was motivated in Q in terms of 
an interaction between the particles of the model, li p > 
1 then the density after the rearrangment tends to be 
distributed equally between sites i and i + 1; if /i < 1 
then the density is more likely to accumulate on just one 
of the sites. The coefhcient A = [R^-'^f'T{2p)/T{pf] is 
determined by the requirement that 

/ dpi / dpi+i Pp(/9i, pi+i |i?) = 1, (3) 
Jo Jo 

which means that all facilitated bonds rearrange with 
unit rate. [Here, T{p) is the usual Gamma function.] 

These dynamical rules respect detailed balance with 
respect to a steady state distribution Psta.t,BBLi{pi}) that 
factorises between sites. In the grand canonical ensemble 
we have 

Pstat,BBL({pJ) ^llPs{p^), Ps{p) CX p^'-^e^", (4) 

i 

normalised so that J^dpPs{p) — 1. Our notation differs 
from 5] in that we use 7 for the Lagrange multiplier con- 
jugate to density, reserving /3 for the inverse temperature 
of the FA model. 

It is clear from Eq. ^ that motion is only possible 
across bonds with fi — 1. We refer to these as facili- 
tated bonds. The steady state contains a finite fraction 
of facilitated bonds, which we denote by 

V ^ in)- (5) 

We also define the mean density, 

P^ip.)- (6) 



Facilitated bonds in the BBL model are the fundamen- 
tal mobility excitations in the system. The interesting 
scaling limit is the one of maximal mean density, p — > 2, 
where facilitated bonds are rare (77 ^ 1). In this limit, 7 
is large, and we have 

p = 2-7-^ + 0(7-'), (7) 
rj = 27exp(-27)^;|^[l + 0(7^')] (8) 

consistent with Q. It was further observed in 0] that 
the dynamics of these excitations in the BBL model can 
be represented schematically by the processes 

01 ^ 11 ^ 10 (9) 

where a 1 represents a facilitated bond (/,; = 1 or fi^i = 
1, respectively), and a an unfacilitated bond {fi — 
or = 0). The above two-step process then produces 
effective diffusion of excitations. When excitations meet, 
they can coagulate via e.g. 101 — > 111 — > Oil — > 010; 
running through the steps in reverse, a single excitation 
can also branch into two. Excitations can never be cre- 
ated unless there is already an excitation present on a 
neighbouring bond, and this is the key motivation for 
the effective FA models presented below. 

When excitations are rare, the rate-limiting step in 
the effective diffusion is the creation of a new excitation, 
01 — > 11. To obtain the typical rate for this process, 
consider a density rearrangement event across bond i + 

Pt,Pt+i,Pt+2 ^ P<,p',+i,p'i+2- (10) 

The process 01 — » 11 occurs when the second bond is 
facilitated in both initial and final states, while the first 
bond is facilitated only in the final state. That is, 

p^+i + Pj+2 = p't+i + p\+i < 2, 
Pi + Pi+i > 2, 
p, + p^+i < 2. (11) 

To work out the typical rate with which these processes 
occur, we should perform a steady-state average over all 
initial configurations with the prescribed mobility config- 
uration (/; = 0, /i+i = 1), corresponding to the first two 
conditions listed. In addition, however, we condition on 
Pi, which strongly infiuences the rate if it is close to 2: 
the third condition given above can then only be met if 
Pijf_i is very small. Thus, we consider the average rate 
for the process (p^ = 2 - ei,f^ = 0, /j+i = 1) ^ {pi = 
— fi — 1, fi+i = 1)1 which can only occur via a den- 
sity rearrangement across bond i -|- 1 as written above. 
The steady-state distribution factorises between sites and 
so we have for this rate, denoted by ri{pi): 
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jldRj^%,+ l Ps(p.+l)Ps(i? - P^+l) 



where we have introduced R — Pi+i + pi+2- The integral 
over p[j^2 ill numerator gives A[p[^^{R — p^+i)]'^^^, 
and the one over p[j^i then a normahzed incomplete Beta 
function B{fi, p.; ei/ R)/ B{fj,, fj,). The remaining average 
over R (and pi+i) becomes concentrated around R — 2 
for large 7, so that in this limit 



i3(/x,M;e./2) _ j;'^''dvvf^-Hl-vr-^ 



ri{ei) 



B{p,p) 



(13) 



Recalling that the local density is = 2 — e^, we note 
that dense sites (those with small €i) lead to small rates 

The relaxation of the BBL model on long time scales 
is determined by sites with small r;. For this reason, 
it is convenient to deduce the distribution of this rate 
from that of pi, or equivalently e^. From Eq. ^ one 
sees for large 7 that the variation of the power law fac- 
tor p^^^ near p — 2 can be neglected, so that Ps(ej) = 
7exp(— 7ei). The typical values of are therefore small, 
£j ~ 7~^, and we can expand the rate as 



(14) 



fi2f'T{p)^ 

Transforming then from the distribution of to gives 



Psie^) 
1 



dri{ei) 



de. 



where 



fin 



(r/r!o 



07 



|(l/M)-lp-(r/ao)l''^ 
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(15) 



(16) 



is a microscopic rate, which acts as an upper cutoff on 
the distribution of rates. [The notation u ^ v means that 
u and V are proportional to each other in the relevant 
limit (large 7).] It is important to note that the small-r 
scaling of the distribution Ps,r(?') ~ r'^'^/f^)-'^ ^ which we 
derived above in the limit 7 — > 00, also applies at finite 
7. This is because the rate for small e; always scales 
as in Eq. pi)) , and the probability density Ps(ei) of 
approaches a constant for small e-i for all 7. 

The long-time behaviour of the BBL model is now con- 
trolled by the behaviour of Ps^r{f) at small r, and, in 
particular, by the exponent p. The time for a mobil- 
ity excitation to diffuse across bond i is of order l/r^. 
The average diffusion time (l/r), with the average taken 
over the distribution Ps,r('")i then shows a change of be- 
haviour at /J, = 1: for /i < 1 it is finite, while for p > 1 



it diverges. This motivates why subdiffusion occurs in 
the second case: for arbitrarily long times t there are a 
significant number of barriers to mobility diffusion that 
have transmission rate <\/t. 



B. Effective FA model 

We now describe the effective model that captures the 
dynamics of the mobile bonds on large length and time 
scales. In this model, the bonds of the BBL model are 
represented by a chain of binary variables {rii}, where 
ni = 1 if the bond between sites i and i -I- 1 of the BBL 
model is mobile, and = otherwise. The variable rii 
corresponds to the BBL variable fi. The process of © 
is then 



(n^ = 0, n,j4 



1) 



1) 



(17) 



It occurs with a rate, rie~^, and the reverse process oc- 
curs with rate r^. Here, determines the concentra- 
tion of sites with Ui = 1, while is a site-dependent rate 
whose fluctuations capture the effect of the fluctuating 
density pi in the BBL model. In our effective model we 
use the convention < < 1; taking a maximal rate 
of unity sets the unit of time. To mimic the distribution 
of rates in the BBL model, we define the disordered FA 
model so that the are distributed independently in the 
steady state, with 



Prin) = il/p)r. 



Ti < 1 



(18) 



in accordance with ([TS]) . 

The rate r.i in the disordered FA model reflects the lo- 
cal density in the BBL model: it is a fluctuating variable. 
In the dynamics of the disordered FA model, we account 
for this fact by randomising when the process corre- 
sponding to (|10p occurs. Hence, we deflne our disordered 
FA model by the dynamical rules: 

rate n,r,e'5("'+^~i'Pann(rO 
{n^,ni+i,n) (1 - rtj,ni+i,r-), 

rate n,+ir,e'^("'^i)Pan„(r-), 



where 



Pi,im{r) (X rPr(r), 



(19) 



(exphcitly, Pann('') = [(1/m) + l]?"^/^), and the vari- 
ables Ui € {0, 1} and < < 1 reside on the sites 
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and bonds of the FA lattice respectively. We ide ntify 
this model as a disordered variant of the FA model [13] , 
since the case Pr{r) = S{r — 1) is the one-spin facilitated 
one-dimensional FA model. We refer to it as the bond- 
disordered FA model since the rates are associated with 
the bonds of the (FA) lattice. 

The 'annealed' distribution of rates after randomisa- 
tion, Pann(?'i), is constructcd such that the model obeys 
detailed balance with respect to 



ri e 



-f3ni 



(20) 



The stationary density of sites with rii = 1 is 

c={n,)^il + e^)-\ (21) 

which plays the part of the parameter rj defined in 
The only other parameter in the model is /u, which cor- 
responds directly with in the BBL model. 

Several other comments are in order. We chose < r < 
1 above, for convenience. As a result, we do not expect 
direct correspondence between time units in the original 
and effective models [comparing (ITSl) and (fTS)) . we have 
effectively set the prefactor to unity] . There is also no 
exact correspondence between the steady states: in the 
FA model, there are no spatial correlations at all between 
the Tii, whereas in the original BBL model neighbouring 
excitations (fi, fi^i) are correlated via the density vari- 
able Pi+i- Finally, also the way rates are linked to 
creation and destruction of excitations does not match 
exactly. In the original BBL model, we saw above that 
the process {f, = OJi+i = 1) ^ {fi = = 1) 

is controlled by the density pi, and is slow when pi 
is close to 2. Translating to the dual lattice of the 
FA model, this corresponds to the controlling rate for 
(jii = 0,ni+i = 1) ^ (rii = = 1) being associated 

with the bond between rti_i and n,;, noi with the bond be- 
tween Hi and TT-i+i as we have posited. Thus e.g. the tran- 
sient appearance of an excitation, 01 11 — » 01, ran- 
domizes ri in our FA model but does not change pi in the 
BBL model so that remains unchanged as well. On the 
other hand, in an effective diffusion step 01 — > 11 — > 10 in 
the BBL model, the second step involves a rearrangement 
across bond i and so a randomisation of pi and hence . 
This is correctly captured in the FA model, and as ef- 
fective diffusion is the key process in the dynamics we 
expect our model to give a qualitatively correct descrip- 
tion of the BBL dynamics. 



C. Model variants 

1. Grand canonical BBL model 

The grand canonical expression (U motivates us to de- 
fine a modified BBL model in which volume is not con- 
served. We use the same dynamical rule (HJ , but replace 



Pp{pi, Pi+i\R) by 

P',{p^,p.+i) = A'(p,p,+i)^-ie'^(''»+'''+^)e(2 - - p,+i) 

(22) 

where the final state is now independent of the volume 
in the initial state. These dynamical rules preserve the 
same equilibrium distribution as that of the original BBL 
model, as given in ([4]). The constant of proportionality A' 
is set by /p dp dp' Pp{p, p') — 1 so that bonds rearrange 
with unit rate, as in the original model. 

We will find that propagation of mobile bonds is sim- 
ilar in models with and without conserved density, al- 
though the relaxation of density fluctuations will clearly 
be different. 



2. Site-disordered FA models 

We also define a site-disordered FA model, in which we 
associate random rates with the sites of the FA chain, 
instead of the bonds. The dynamical rules arc then 



(n,, 1 - nj+i,r,,r-+i), 

rate n,r,+ie'^("'+i-i)Pann(r^+i) 
(1 - ni,ni+i,r[,ri+i), 

rate n,+me>"^''^-'^ P.nnir',) 



As for the bond-disordered FA model, also this model 
does not exactly capture how rates are linked to rear- 
rangements in the BBL model; but it does provide for 
rates to be randomised every time an excitation makes 
an effective diffusion step, which is the key property for 
the physics. Indeed, we will see below that the excitations 
behave similarly for bond and site disorder. However, on 
introducing probe particles to these disordered FA mod- 
els, one finds that the site-disordered model provides a 
better match to the BBL model dynamics. The reasons 
for this will be explained below. 



3. FA models with quenched disorder 

Finally, it is convenient to define FA models with 
quenched disorder, in which the rates r^ do not depend 
on time. The distribution of rates is simply Pr{ri) in 
that case. Interestingly, we will find that quenching the 
disorder in this way has very little effect on dynamical 
correlations (after disorder averaging). Wc note that the 
quenched bond-disordered FA model has a mapping to 
a disordered model of appearing and annihilating defects 
( AA model) , and inherits from the latter an exact duality 
mapping, as in the pure case [21 1. 



III. MOBILITY EXCITATIONS 

We now consider the dynamics of the mobility excita- 
tions in the BBL model, always in the interesting limit 



5 



1 

0.8 
0.6 
0.4 
0.2 






i 




■ 

1 


4 












V 


[ / 




■ I 






\ t 




1 i 





200 400 600 800 1000 
site 



1 

0.8 
0.6 
0.4 
0.2 




200 400 600 800 1000 
Site 




200 400 600 800 1000 
site 



FIG. 1: Typical realisations of the propensity, with times such that the spatially averaged persistence function satisfies 1 — 
{p{t)) ~ 0.4. Large values of the propensity indicate sites that are very likely to have relaxed, on this time scale. The models 
with subdiffusive dynamics have large jumps in the propensity, which arise from sites that relax very infrequently. (Left) 
BBL model at /i = 2 and -q = 0.01. (Center) Bond-disordered FA model ai fi = 2 and c — 0.01. (Right) Pure FA model at 
c — 0.01. In the pure case, there are no barriers for excitation propagation, and the propensity is smooth [except near sites 
with n,;(t = 0) = 1, where it is maximal]. 



where p is close to two. First consider the regime where 
the parameter fi is small (much less than unity). The 
BBL model in its steady state then has a bimodal distri- 
bution of densities with sharp peaks near zero and two. 
In that case, it describes diffusing vacancies in a one- 
dimensional soUd (i.e. high-density background). The 
disordered FA models, on the other hand, all reduce to 
the pure FA model in the limit of small fi. All models 
then exhibit dynamical scaling when the excitation den- 
sity is small, with exponents 

(z,;/) = (2,1), fi<l (23) 

Here z is the dynamical exponent that sets the relative 
scaling of space and time, while the correlation length 
scales as the average distance between excitations, i.e. 
5 ~ ry""^ for the BBL model and ^ ^ for the FA case. 

However, the case of /i > 1 is qualitatively different 
from that of small /i. For example, the mean time as- 
sociated with rearrangements is (1/r) which diverges for 

> 1 as explained above [recall equations ([15]) and p8|) ]. 
We therefore expect the disorder to have a strong effect: 
this is clear from plots of the propensity [2^, which we 
define in terms of the persistence function, Pi{t). This 
function takes a value of unity if the state of site i has 
not changed between time zero and time t, and Pi{t) = 
otherwise. The (time dependent) propensity for a given 
initial condition of the system is then 1 — {pi{t)) dyn, where 
the average is over the stochastic dynamics of the system, 
but with the initial condition fixed [l^l ■ We show sample 
plots in figure [TJ In the BBL model, sites with density 
close to two act as barriers to propagation of mobility; 
in the FA model the same effect arises from bonds with 
small rates. 

A more quantitative measure of the effect of the disor- 
der is its effect on dynamical scaling. It was observed in 
that the BBL model has scaling exponents close to 

(z,z.) = (l + /i,l), n>l. (24) 

An effective model of a single excitation propagating in 
a quenched environment of random energy barriers gives 



this scaling, if the distribution of rates for crossing the 
barriers is Pr{r) [2^. We will discuss below why this 
quenched result is applicable to the BBL model, which 
has no quenched disorder. First, though, we show that 
this subdiffusive scaling can be observed in the four-point 
functions of both BBL and disordered FA models. 
We consider the correlation function 

Gi{x,t) = {6p,+,{t)6pM) (25) 

where Spi{t) = Pi{t) — {pi{t)); Pi{t) is the persistence oper- 
ator defined above, and the averages are over both initial 
conditions and the stochastic dynamics. The normalised 
four point susceptibility is 

X4„W =5]G4(a;,i)/G4(0,0 (26) 

X 

In one dimension X4n(i) is a direct measurement of a 
growing length scale, when normalised in this way (20l . 
[24| . (Note that we evaluate averages in an ensemble with 
fixed 'chemical potential' 7, so that the mean density p is 
allowed to fluctuate. Since the dynamics conserve p, this 
choice does affect the value of X4n(i) [I]-) In the scaling 
limit (p — > 2 from below) , we expect 

X4„(i) ^ {Qty^^fint/o (27) 

where ^ is the correlation length whose scaling was given 
above, z is given by ([23|l or ([24|l . as appropriate, ft is 
a microscopic rate and f{x) is a scaling function that is 
constant at small x and decreases as x~^/^ for large x. 
We argue below that the for the BBL model, the rate 
is equal to fio [recall (fT5|) ]. while for the FA model, we 
have ~ c (for small excitation density c). 

We show results in figure [2l Both models are consis- 
tent with ((27|l : we also find that the FA model exhibits 
the same scaling as the BBL model, if we identify the 
excitation densities c and rj. Hence we argue that the 
disordered FA model is an appropriate effective theory 
for the BBL model. 

Figure [2] demonstrates further that BBL models with 
and without conserved density behave very similarly. We 
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FIG. 2: Subdiffusion in the FA and BBL models is appar- 
ent in the four point susceptibility. (Top) We show X4n(t) in 
the BBL model at various ^ and 77. For the case ^ = Q, 
we use a binary distribution of on-site densities Pb{p) oc 
aS{p) + 5{2 - p - 0+), so that r; = a(2 + a) /{I + af; the 
small offset 0^ ensures that the facilitation constraint is well 
defined. The dashed lines show the power law predictions of 
(|24[) . and simple diffusive scaling (|23|) for fi — 0. At large 
times, X4n(i) saturates at a value of order r]~^: this is appar- 
ent in the data for p = 0. The label (G) denotes data for the 
grand canonical variant of the BBL model. (Bottom) Bond- 
disordered FA model data with the same values of p, showing 
that this effective model captures the four point correlations 
of the BBL model. The labels (Q) and (S) denote data from 
the model with quenched bond disorder and fluctuating site 
disorder respectively; the different variants have very similar 
behaviour. 

conclude that the conservation of density is not relevant 
for scaling: this is consistent with the use of the disor- 
dered FA model as an effective theory, since that model 
has no conserved density. Quenching the disorder in 
the FA model has only very weak effects on disorder- 
averaged properties such as X4n(i); finally, differences 
between bond-disordered and site-disordered models are 
also very small. 



A. Effective barrier models for single excitations 

Figure [2] is clear evidence that both BBL and disor- 
dered FA models have excitations that propagate subd- 
iffusively at large /i. Further, the dynamical exponents 
for all the models with p, > 1 seem to satisfy (1^^ . 



For the FA model with quenched bond disorder, this 
result is to be expected since the motion of indepen- 
dent random walkers in this kind of environment is well- 
understood [2^ and does indeed satisfy ([M]) . However, 
it was argued in that fluctuating disorder should lead 
to 

z = fi, p>2 (28) 

and z — 2 otherwise. This result is inconsistent with the 
data. 

We are not aware of any analysis of fluctuating disor- 
der that is slaved to the motion of the random walker. 
In this section, we give an argument that explains the 
applicability of ([M)) to the FA model with fluctuating 
disorder, and hence to the BBL model. While this is not 
a rigorous proof, the various stages of the argument have 
been verified by direct simulation. 

To describe the motion of a single excitation in a 
disordered environment, we consider a simple barrier 
model [23]. A single particle moves on a chain of 
sites, with independent random hop rates {ri} on the 
bonds, distributed according to Pr{r)- We consider both 
quenched and fluctuating disorder: if the disorder is fluc- 
tuating, then each random rate is redrawn from the dis- 
tribution Pannl?") whcn the bond is traversed by the ran- 
dom walker. We have verified by simulation that both 
variants of this model do indeed satisfy We ex- 

plain this result using an argument related to that of 
le Doussal, Monthus and Fisher 25]. The effective dy- 
namics scheme that we use for the barrier model with 
quenched disorder was described in ^26], where it was 
shown that the effective dynamics are a good description 
of the quenched barrier model, as long as the exponent 
p, is large. We give a brief description of the effective 
dynamics here, referring to [26j for details. 

In ^2^] , the authors proposed an effective dynamics for 
a random walker in a (quenched) one-dimensional energy 
landscape, made up of 'barriers' and 'valleys'. At each 
stage of the effective dynamics, the smallest barrier in the 
system is removed, and the particle moves to the bottom 
of the valley that contains the origin. The time associ- 
ated with this process is the inverse transmission rate of 
the barrier that was removed. For models in which the 
energy landscape has short-ranged correlations, this ef- 
fective dynamics mimics the real dynamics of the random 
walker. 

For the quenched barrier model, every site is at zero 
energy, and they are separated by barriers of varying 
heights. The effective dynamics involves successive re- 
moval of the smallest barriers. Thus, at a given stage of 
the dynamics, the remaining barriers divide the system 
into 'effective traps'. As discussed in [2^, the barrier 
model requires a modification to the scheme of [25| . in 
that the time at which barrier i is removed depends both 
on the rate and on the widths of the effective traps 
to the left and right of barrier i. If the widths of these 
traps are h and I2, the time associated with barrier i 
is determined by t^^ = ri{l/li + l/h)- 
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To arrive at the subdifFusive scaUng of the quenched 
barrier model, we assume that, at time t, effective traps 
have a typical width i{t). All barriers with t.; < t 
have been removed; typically, these barriers have ^ 
£{t)/t [28j]. Thus, the density of remaining barriers is 



rm/t 

£{ty^ ~ / drPr{r) 
Jo 

which yields (for large t) 



(29) 



(30) 



The mean square displacement of the diffusing excitation 
scales with ^(i)^, so we identify the dynamic exponent 
z = 1 + n, consistent with ([M)) . Recall that the effective 
dynamics scheme applies only for large fi, so we do not 
recover the diffusive result, z = 2 for ^ < 1. 

We now apply this scheme to models with fluctuating 
disorder. The fluctuations in the disorder have two main 
effects. Firstly, once large barriers have been crossed, 
their rates are randomised. Thus, if multiple crossings 
of the same large barrier are important for the quenched 
model, we expect different behaviour for fluctuating dis- 
order. However, a central assumption of the effective dy- 
namics is that the time to travel a distance £ is dominated 
by the time required for the first crossing of the largest 
barrier between the initial and flnal sites [2^. Hence, 
multiple crossings of large barriers are ignored in the ef- 
fective dynamics, which should therefore be consistent 
with fluctuating disorder. The second effect of fluctuat- 
ing disorder is that barrier transmission rates are being 
randomised as the excitation moves around, so a barrier 
which previously had a large transmission rate may ac- 
quire a new rate that is very small. This new rate would 
then act as a high barrier and so have a strong effect on 
the resulting motion. As before, the time taken to move 
a distance £ will be given by the time taken to cross the 
largest barrier between initial and flnal states: this might 
be a barrier that was present initially, or one that ap- 
peared as the excitation moved through the system. The 
key point here is that the system is in a steady state, 
so the introduction of new barriers occurs with the same 
rate as the removal of barriers of the same size. Since 
barriers are removed only when they are crossed, barriers 
that appear in the system are typically of a size compa- 
rable with those that have already been crossed at least 
once. In the language of the effective dynamics, these 
barriers are 'irrelevant'. For these reasons, the effective 
dynamics apply equally well to models with quenched 
and fluctuating disorder, and we expect the dynamical 
exponent z = 1 + ii for both cases. 

Of course, the situation would be very different if the 
fluctuating disorder was annealed in a two-sided way, 
where every time an excitation moved to a new site one 
randomises the rates for both of the barriers adjacent to 
that site. This would produce a continuous-time random 
walk 23], with subdiffusion exponents as in Eq. 



To obtain the scaling of length and time scales in the 
FA and BBL models, we note that the equilibrium spac- 
ing between defects sets the dynamical correlation length 
^ (for these one-dimensional models). We deflne the per- 
sistence time Tp by (pi(Tp)) = 1/e, where Pi{t) is the 
persistence function, deflned above. As in the pure FA 
model, Tp scales with the time taken for an excitation 
to propagate a distance ^, so we identify Tp ~ 
consistent with (|27|) . To obtain the scaling of with 
the excitation density rj or c, it is useful to rephrase the 
scaling argument associated with the effective dynamics. 
In the BBL model, Eq. implies that the fraction of 
sites with rate r, < r scales as 



n(r) ~ {r/noY/^ 



(31) 



for small r. Thus, moving a distance £ typically requires 
the particle to cross a barrier whose transmission rate is 
Vi ~ no£~^. The time taken to cross such a barrier is 
'i£i+tj.^ This is consistent with 







typically n ^ £/r, 

and it allows us to identify the coefflcient in P?)) 
with JIq in In the FA model, crossing a barrier 

with transmission rate typically requires a spin to flip 
from state to state 1, and this process occurs with rate 
~ c. Thus, moving a distance £ typically requires the 
crossing of a barrier with ~ £~'^, which takes a time 
Ti ~ £/{ric) ^ c~^£^^^^. Thus, we identify the coefficient 
ri in (|27p with the inverse excitation density c. Overall, 
for /X > 1, we arrive at Tp ^ c~'^~^ for the FA model, 
and Tp ~ \h\{\/lrj)\~^^r]~^~^^ for the BBL model, where 
we used flo ~ ~ [ln(l/f7)]~^- 

We conclude that the effective dynamics scheme pre- 
sented here captures the propagation of mobility exci- 
tations in the BBL and disordered FA models on large 
length and time scales, even though the BBL model 
has non-trivial dynamical correlations in the densities pi 
which the coarse-grained FA model neglects. This analy- 
sis demonstrates that the scaling properties of the persis- 
tence time and the four-point susceptibility, as expressed 
in and ([?f|) . can be understood in terms of in- 

dependently propagating (non-cooperative) excitations. 



B. Long-time limit 



So far, we have considered time scales up to the per- 
sistence time Tpi excitations move distances smaller than 
their typical spacing, and can be treated independently. 
We now turn to much longer time scales. The assump- 
tion of independently propagating defects in one dimen- 
sion leads to a persistence function consistent with the 
results of Q: 



p{t) = {p,{t)) = exp 



(nt) 



l/z 



(32) 



for i > 1. 

However, for times larger than Tp, this prediction fails. 
For example, in the site-disordered FA model, the frac- 
tion of sites with rate < r is n{r) = r^/^. At infinite 
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temperature (/3 = 0, c = 1/2), the facilitation constraint 
in the FA model can be ignored (most sites are faciliated). 
In that case, the typical time taken to flip (for the first 
time) a site with initial rate is = l/r^. Thus, the per- 
sistence function decays as {pi{t))p=o ~ n{t~^) ~ t~^/^. 
Lowering the temperature in the FA model only slows 
down the dynamics, so {pi{t)) ^ t~^/>^ for all times and 
temperatures. Thus, (j32p must break down at long times: 
we attribute this breakdown to the fact a single site with 
a small rate can block the motion of several excitations. 

We now consider this long-time regime in more de- 
tail, and return to the effective dynamics picture, work- 
ing with a finite density of excitations, 77. If the density 
of relevant barriers is larger than the density of excita- 
tions, each effective trap typically contains at most one 
excitation, and excitations can be treated independently. 
However, when the spacing between relevant barriers be- 
comes larger than the distance between excitations, one 
enters a different regime. To see this, note that the 
typical time scale associated with rearrangement of a 
'slow' (relevant) site i in the BBL model is generically 
Ti = \/[ri{rii-i -|-77i-|_i)] where ?7i_i and rji+i are the exci- 
tation densities in the effective traps to left and right of 
site i. (Recall that is the rate with which the relevant 
site rearranges, given that there is an excitation adja- 
cent to that site. Thus, the time taken to flip a relevant 
site depends on the density of excitations in the adjacent 
traps.) 

In the short-time regime where there are many more 
effective traps than there are excitations, then we can 
write ?7i ~ 1 /h if trap i contains an excitation, and 77^ = 
otherwise (as above, h is the width of effective trap i). 
Considering a site i for which one of the adjacent traps 
contains an excitation, we arrive at the scaling relation 
Ti ~ h/ri ~ (-/rt, as discussed above. However, if there 
are more excitations than effective traps, we expect the 
density in each trap to be close to its equilibrium value 
Tji ~ 77j_|_i ~ 77. Thus, we expect ^ {2rjri)^^ for £77 ^ 1. 
In both cases, for a given time we use 7i(r) ~ (r/57)^/^ 
to evaluate the fraction of sites with Ti < t: the mean 
spacing between these "relevant" sites is i{t). The result 
is 



(^7^/e)l/^ 



m « e 
m » ^ 



(33) 



where we have used C ~ '? ^- For long times, the ex- 
ponent i sets the time dependence of £{t): this result 
applies for all fj, > 0. In the short time regime, (|33p is 
consistent with the analysis of the previous section, and 
with Eq. (P7| . as long as /i > 1. However, if < 1, the 
motion of excitations is diffusive: thus, in the short time 
regime, there is no distinction between relevant and ir- 
relevant barriers. This means that if /i < 1, the spacing 
between relevant barriers, i{t), is only well-defined in the 
long time limit, and the short-time scaling regime of ((33|) 
does not exist. 

We observe that in the long time limit, £{t) represents 
the mean spacing between isolated sites with small rates 
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FIG. 3: (Top) Persistence function, plotted against the scal- 
ing variable At in the site-disordered FA model. We show 
data for that c = 0.5, 0.27, 0.12, 0.047, 0.018 (decreasing from 
left to right, these are inverse temperatures /9 = 0, 1, 2, 3, 4). 
As discussed in the text, the persistence time scales as Tp ~ 
c~^~'^ ~ c~*, but we use the variable variable Xt to verify the 
long time prediction of (|35|) : the dashe lines shows a power 
law with exponent i. (Bottom) Persistence at infinite tem- 
perature and varying /i. Again, dashed lines show power laws 
with exponents predicted by psp . 



Tj, and these sites dominate the long-time limit of the 
persistence function. That is, in the long time scaling 
regime, ([5^ is replaced by 



p{t) £{t)-' c (m/0-'/^ 



(34) 



The crossover between the two scaling regimes occurs 
when £(t) ~ ^ ~ t;"^ or c~^, respectively. This can be 
observed in the long-time behaviour of the persistence 
function in the site-disordered FA model. To obtain the 
long-time limit of this function more quantitatively, we 
decompose the persistence p{t) = cpi{t) + (1 — c)po{t) 
into contributions from sites that were initially in states 
1 and [these two populations have weights c and (1 — c) 
respectively]. Then, in the long time regime ^ ^ ^, we 
have Ti ~ (2cri)~^ for sites that were initially in state 
1, given that each of the two facilitating neighbour sites 
contains a defect with probability c. For those sites ini- 
tially in state 0, the spin flip rate is suppressed by , 
leading to Ti ~ (2ce~''ri)^^. The persistence functions 
are then estimated as the density of sites with Ti > t, 
giving pi(t) = 7i(l/(2ct)) andpo(i) = 7i(l/(2ce-''<)), re- 
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spectively. Using n(r) = r^^^ , we thus arrive at 

p(t)~{\t)-^/^' ^(t)»e~c-i (35) 

with 

A = 2ce-^[l + c(e-'3/^-l)]-^ (36) 

In the hmit of dilute excitations (c <C 1) this reduces 
to A ~ c^. Identifying p{t) with l{t)^^ , this resuh is 
consistent with ((33l) . since we argued in Section IIII Al 
that c for the FA modeL 

Results are shown in Fig.[3l at infinite temperature the 
power-law behaviour of the persistence is clear. At lower 
temperatures, the crossover to power-law behaviour oc- 
curs deep in the tails of the persistence [p{t) ^ c]. In the 
FA model, this long time regime can be demonstrated 
by simulations at high temperature, on relatively short 
time scales. However, in the BBL model, the equivalent 
of the high-temperature regime requires small 7, increas- 
ing the prefactor f2oj ^^'^ reducing the fraction of sites 
with small r^. This means that very long simulations are 
required to access the long-time limit in the BBL model, 
and we do not show numerical data in this case. However, 
the simulations of the site-disordered FA model confirm 
the validity of the arguments of this section, which apply 
to both FA and BBL models. (To observe the long-time 
regime in the bond-disordered FA model, one would need 
to define and measure a persistence observable on bond i, 
associated with the rearrangement of density across that 
bond.) 

IV. PROBE PARTICLES 

It is a familiar feature of kinetically constrained models 
that propagation of probe particles is different from that 
of excitations [3, 4, 6, 27]. We now turn to probe particle 
motion in the BBL and disordered FA models. 



A. Probes in the BBL model 

We introduce (non-interacting) probe particles to the 
BBL model as follows. A probe can move along a bond 
when density rearranges across that bond. If the bond 
connects sites i and i -I- 1, then after the rearrangement, 
the probe occupies site i with probability pi/ {pi -\- Pi+i)- 
The joint stationary distribution for the probe position 
and the BBL densities is 

-Pprobe,stat(-^, {Pi}) CX Pi^JsT J]^ -Ps (Pi ) ) (37) 

i 

where X is the position (site index) of the probe. Thus, 
the probability of finding a probe on site j = X is is pro- 
portional to the local BBL density on that site, pi=x- 
This is consistent with the probe representing a typical 
particle in the BBL model, before the coarse-graining 



into the densities pi is carried out. An alternative rule, 
which is more consistent with the effective FA model de- 
scribed below, is to assign a probe with equal probability 
to sites i and i + 1. As we discuss below, the excita- 
tion motion in these models sets bounds on the motion 
of the probes: we are primarily concerned with the situ- 
ation in which these bounds are saturated, in which case 
details of the microscopic probe motion should be irrel- 
evant. When the bounds are not saturated, the choice 
of dynamical rule does produce quantitative differences, 
although qualitative features are preserved. 



B. Probes in the FA model 

We couple probe particles to the FA model using the 
method of Probes can hop between pairs of adjacent 
sites only when both sites have n = 1; they attempt these 
hops with unit rate. With these rules, the equilibrium 
distribution analogous to ([57]) is independent of X: that 
is, the distribution of the probe position decouples from 
the excitation variables rii and the rates r^. 

From the data presented above on the excitation dy- 
namics of the site-disordered and bond-disordered FA 
models, one might expect that the two model variants 
also exhibit similar probe dynamics. However, this is not 
the case because barriers to excitation diffusion act dif- 
ferently on the probes. To see this, consider the site- 
disordered model, and suppose that site i starts with 
rii = and with a small rate r^. The probe cannot cross 
this site until its excitation state changes to = 1. The 
rate for this is of order cr^, and so the rate for a probe 
to cross this site also vanishes with r^: high barriers for 
excitations (small r^) are also high barriers for probes in 
the site-disordered FA model. 

Now consider the bond-disordered FA model, focusing 
on a particular bond i, with a small rate r^. The probe 
particle can cross this bond ii rii — \ and n^+i = 1: this 
state can occur on time scales much shorter than (cr^)"^ 
if an excitation arriving from the right facilitates n^+i 
and another excitation arriving from the left facilitates 
Ui. This process sets a rate for crossing the slow bond 
that is independent of rj. So the barriers for excitation 
diffusion have a much smaller effect on probe propagation 
in the bond-disordered model. (One way to avoid this 
behaviour would be to allow the probe to move along 
bond i only when the rate for that bond is randomised, 
but we have not pursued this as we wanted to keep the 
probe dynamics similar to that used in [3].) 

It is clear that in the original BBL model, the barriers 
for excitation diffusion do also act on probes. A high 
barrier here is a site with density pi « 2. This can take 
part in a rearrangement only once a rearrangement of 
neighbouring sites has produced a low-density < 
2 — Pi or yOi_|_i < 2 — Pi. The rate for these processes, 
and hence for probe diffusion across site i, vanishes as 
Pi — > 2. In summary, only the site-disordered FA model 
can provide an accurate representation of the BBL probe 
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FIG. 4: Comparison of root mean square probe displacement 
^rms(t) = {Xitf')^!'^ with the inverse of the persistence func- 
tion l/p(t), for the site-disordered FA model. We show data 
at ^ = 3 (top panel) and = 4 (bottom panel). We have 
offset the data at c = 0.25 for clarity: the ordinate is lOt in 
that case. In the long-time limit, the persistence scales as 
t-^l^ (this power law is illustrated with dashed lines). At 
long times the persistence sets a bound on the scaling of the 
mean square displacement, which appears to saturate for long 
times. 



dynamics because it retains the effect of high barriers 
on the probes. We therefore do not consider the bond- 
disordered case in the following. 



C. Results for probe motion 

In the BBL and site-disordered FA models, the preced- 
ing discussion illustrates that sites with small rate are 
able to block the propagation of probes. Taking the FA 
model for concreteness, the probe cannot pass any site 
for which = for all times between and t. As dis- 
cussed in Section fill Bl the mean spacing between these 
sites scales as t^/^ for large times t. This sets a limit on 
probe motion 



FIG. 5: Mean-square probe displacements in the BBL and 
FA models, for times shorter than Tp. (Top) BBL model. 



(38) 



For ^ < 2, this bound is irrelevant: the probes simply 
diffuse. However, for /i > 2, we expect this bound to 
be saturated at large times. At long times, we have 
e{t) ~ p(i)"\ [recah JM])]. Thus, Fig. H demonstrates 
that the bound ([55)1 does saturate at long times, although 



For fj, > 2, the dashed lines show the predictions of (|40|l : for 
/i = 1, we show a diffusive law {X^{t)) ~ t. (Bottom) Site- 
disordered FA model, with the same power-law predictions. 
In both cases, we expect a crossover at long times to the 
asymptotic scaling of (|38|l . 



we note that the times required are quite large, even at in- 
finite temperature (c = 1/2). Physically, the length scale 
£{t) represents the size of an effective trap: saturation of 
the bound requires that the probe particle explores the 
whole of the trap before the barriers delimiting the trap 
become irrelevant. The scaling arguments presented here 
do not allow us to estimate the time required to reach 
this regime. However, Eq. ([38|l shows that probe propa- 
gation must be asympotically subdiffusive for all /i > 2, 
and the data are consistent with saturation of this bound 
throughout this regime. 

We emphasise that while Fig. [4] demonstrates that ([38]) 
holds on long time scales in the FA model, the scaling 
arguments presented here apply equally well to the BBL 
model, so asymptotic probe motion in that model must 
be subdiffusive for /i > 2. 

We now turn to time scales shorter than Tp, for which 
the length scale £ again sets a bound on probe motion. 
Following we decompose the probes into two popu- 
lations: those that have moved at least once, and those 
which have not moved at all. We denote the fraction of 
probes that have moved at least once by 7r(i). Confine- 
ment of probes by sites that are persistently unfaciliated 
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again sets an upper bound on the displacement of probes: 



1 - 7T{t) 



< 



i{ty 



(39) 



where the left hand side is the nth moment of the dis- 
tance moved by those probes that have moved at least 
once, and the scaling of £{t) was given in p3|) . Again, 
saturation of this bound occurs when motion of the probe 
particle within the effective trap is fast enough that the 
the probe can delocalise within the trap before the bar- 
riers delimiting the trap become irrelevant. In the joint 
limit of large time and large /i, the time scales associated 
with adjacent barriers become well-separated al- 
lowing equilibration to take place. Thus, for large fj,, we 
expect the bound of ([39| to be saturated for times t ^ 1, 
even if t <C Tp. 

Assuming saturation of the bound ([5^ and t <C Tp, the 
probe persistence scales as 1 — 7r(<) ~ 1 — exp(— £(t)/^) ~ 
£{t)/(^. (This is the same scaling as for the excitation 
persistence in ([32]).) Combining this with (|33p . we arrive 
at 



(40) 



Our simulations are restricted to finite time scales and 
values of fi that are not very large, so we are not able 
to investigate this bound in detail. However, the results 
shown in Fig. [5] are certainly consistent with the predic- 
tion of (gni). 



V. CONCLUSION 

To summarise our main results, we have established 
that probes and mobility excitations both propagate sub- 
diffusively in the BBL model, and that this subdiffusive 
behaviour can be reproduced in a simple effective FA 
model. This observation allows us to analyse the subd- 
iffusive motion, and to predict the dynamical exponents 
for both excitations and probes in the subdiffusive regime 
[Eqs. ([M]) and ([55)1 ]. A key part of the reasoning consists 
in showing that quenched and annealed disorder lead to 
qualitatively the same behaviour. This allowed us to de- 
duce that correlation length and time scales are related in 
these models, by r ~ When /i is large, we conclude 

that the very broad distribution of rates in these mod- 
els leads to a relaxation time that increases much more 
quickly than the associated length scales, on approaching 
the glass transition. 

We also identify two kinds of subdiffusive motion in 
these models. On time scales 1 ^ i ^ Tp, mobility 
excitations propagate independently and subdiffusively, 
according to 



(X^) ~ (fit) — 



(41) 



One has to remember though that this {X^y^"^ does not 
define a lengthscale for probe motion, since it arises from 



an average over a dominant population of probes that 
have not yet moved, and a smaller population that has 
moved by (ilt)^/(^+^\ For the same reason the expo- 
nent for the scaling of {\X{t)\") in Eq. (|40p is not simply 
proportional to n. On time scales t ^ Tp, excitations 
coagulate and branch, and it is not consistent to discuss 
motion of a single excitation. However, in this long-time 
regime, probe particles propagate subdiffusively, accord- 
ing to 



(42) 



The presence of different dynamical exponents for probes 
and excitations may seem surprising, but we emphasise 
that ([4T|) and ([42]) apply in separate scaling regimes. 
(When the concentration of excitations is small, the per- 
sistence time Tp 1 separates two well-defined scaling 
regimes; of course t is always taken to be large compared 
to unity.) 

Conceptually, it is interesting to note that in the pure 
FA model at low temperature, relaxation is controlled by 
rare active sites (defects). In the disordered model, on 
the other hand, rare inactive regions (sites with small r^) 
play at least as important a role. 

Finally, our results forprobe particles imply that the 
Stokes-Einstein relation [9| between relaxation time and 
probe diffusion constant, Dt ~ 1, has broken down com- 
pletely in these systems. In the pure FA model, Dt di- 
verges at low temperatures On the other hand, in 
the disordered model, the presence of sites (or barriers) 
with arbitrarily small rate means that the persistence 
decays as a power law for large times, while the motion 
of the probes is subdiffusive even in the long time limit. 
However, we can define an analogue of the Fickian length 
£p — (X^(rp))^/^ which represents the distance travelled 
by a probe, through repeated encounters with a single 
excitation [j| . If the bound of ([59)1 is saturated we arrive 
at ip ~ £{Tp) ~ For the site-disordered FA model, this 
leads to £f ~ c~^, at least for large ^; on the other hand, 
in the pure FA model, ^ c~^/^. Physically, confine- 
ment of the excitation in an effective trap means that it 
facilitates any probes in that trap very many times, al- 
lowing the probe to delocalise thoughout the trap. In this 
way, the presence of large barriers to excitation diffusion 
in the BBL and disordered FA models strengthens the ef- 
fects discussed in [1, |3] , in which the square of the Fickian 
length represents the number of hops that a probe makes 
through multiple encounters with a single excitation. 
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